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TECHNICAL MEMORANDUM 7Kir,l 


GRAVITATION 
I. INTRODUCTION 


During the period June <1 to August 19, 1977, investigations were conducted 
by the author in the Cryogenic Physics Branch of the Space Physics Division of 
the Space Sciences Laboratory, Marshall Space Flight Center. This branch is 
involved in preparing tho Stanford gyroscope experiment to detect the gravita- 
tional induction fields generated by rotating bodios l 1|. These investigations 
wen*, therefore, in tho general area of consideration of gravitational fields 
predicted by various currently viable, or proposed, theories of gravitation. 

The question of the existence of black holes was treated with the aid of 
the symbolic manipulation language MACSVMA in collaboration with Dr. Richard 
Pavollo of I.ogicon Ine. , a consultant to Hanscomb Fields's Air Force Cambridge 
Research Laboratories and to the Redstone Arsenal. Use of the term "black 
hole" will be reserved for later discussion, but trapped-surfaces [2), photon 
capture orbits, and photon capture impact parameters (3) were indeed found to 
exist in the gravitation theories of Einstein |3|, Hrans-Dieko | 1 1 , Rosen [5|, 
Lightman-I.ee (0), and Yang |7|. These properties strongly imply the existence 
of black holes in any gravitation theory with a Riemannian geometry. 

A two-tensor theory of Yilmaz l M— 10 ) with scalar field was completely 
and exhaustively studied. It was determined to be inconsistent, noneovariant, 
and completely devoid of gravitational effects. An earlier argument of 
Kraichnan [11) indicates that the result should lie true for a wide class of 
gravitation theories. 

Anisotropic cosmological models were investigated in Rosen's theory 
(12,13) of gravitation. It was found that the simple cosmologies liecome iso- 
tropic much more rapidly than in Kinstoin's theory. The nine Bianchi classifica- 
tions have been extended to Rosen's theory for computations in the spatially 
homogeneous models. 

Inhomogeneous cosmological models were investigated in the Einstein 
theory. The approach was to generalize* spatially homogeneous models by 
raising the dimension of the invariance group by one parameter in the reverse 
of Inonu-W igner's contraction | 1 1| . Three eases are considered; 



T 


a. The Introduction of a peculiar velocity into self-similar cosmologies 
1 15| to sm' what the major changes in observations would 1** (the generalization 
of Bianchi-Typo Mil to the inhomogeneous case). 

t>. The generalization of Class III locally rotationally symmetric 
universes (1G| to inhomogeneous spaces after a method suggested by Wainw right 
117|. 


c. The extension of spatially homogeneous universes to inhomogeneous 
ones by adding an inhomogeneous nonsynch ronous term to the metric. 

The Ilianchi groups operate on observer rest spaces which slice through the 
inhomogeneous space sections. This is a continuation of earlier work by the 
author 1 1«| and is a view of cosmologies which has been discussed in a 
different form by Spero and Baierlein 1 19| . 

A solution has lieen found for cosmological models in Yang's theory of 
gravitation which is not an Kinstein space. The observational consequences are 
straightforward to investigate since the metric is similar to that of the steady- 
state spatially homogeneous cosmological models. 

The Ncwman-Pcnrosc formalism [20] for It iemannir.n geometries allows 
one to investigate space time's characterized by the presence or absence of one 
property or another. The use of the formalism to find solutions to any viable 
gravitation theory is being investigated. The relevant projierty is that the 
gravitational field propagate on twisting null rays. Thus, it is hoped that the 
means may lx* found for obtaining solutions appropriate to the exterior of a 
rotating star and to universes containing circularly polarized gravitational waves. 
Progress in those areas is outlined in this report. 


II. FXISTFNCE OF BLACK HOLES IN 
GRAVITATION THEORIES 


Black holes have lieen an accepted feature of general relativity for 
several years, but only within the last 10 years have they become of widespread 
interest in gravitation and astrophysics [3|. Because they contain a singularity 
of zx'ro volume and infinite matter <lensity and, therefore, represent regions 
from which light cannot escape, they are a pathology of the theory, indeed, 
many relativists have come to regard their existence as an indication of a flaw 
in general relativity. They believe that black holes have become a symptom of 
an intrinsic disease in general relativity itself for which then* is no cum. 
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This situation lias led, in Homo cases, to attempts to find altornativo 
theories of gravity which will not predict the existence of black holes (5, G| . It 
has also sparked the hope Uiat a proper marriage of general relativity with 
quantum field theory would prevent black -hole formation in an evaporative 
mass-loss process (21 1. 'I' he mass-loss rate would conceivably match the 
collapse rate in such a way as to avoid the formation of the black holo. 

This section will show that the same situation exists in all presently 
known viable gravitation theories. The criteria for the existence of black holes 
will be presented and then applied to these gravity theories. Examination of 
the two-dimensional collapse scenario with scalar quantum field theory leads 
fimdly to the conclusion that black-hole eva|x>ration will lie predicted by all 
these theories, implying no distinction between them and general relativity. 

I/?t a gravitation theory define matter trajectories by the geodesic 
equation 


'a lib f 


0 


where II indicates covariant differentiation and t is tangent to the trajectory. 
Then a congruence of trajectories w ill obey the geodesic deviation equation 

relating the connecting vector n J to the Hiemann tensor It' 1 , • 

bed 




t 


bed 
n f 


Let these trajectories lx- cmlieddcd in a Hicmannian manifold or. which the 
Hianchi identities hold. Then the effects of the deviation equation w ill lie 
contained in the N'c wm an- lH? n rose equations |20|. These equations connect the 
evolution of the geometric d spin coefficients to each other. The spin coefficients, 
in turn, contain the information about the physical Ix-havior of the matter 
trajectories' congruences. 

The spin coefficient of interest is p , the complex expansion of the 
congruence, describing the expansion 
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m iH I i ho rotation 

of the rays (the brnekot, | |, denotes antisymmefrl/ation) 
p «= 0 + i:o . 


Lot f l*o tangent to an outgoing null geodesic ami n taneopt to an evoing null 
a a 

geode nic. Let the complex vector n span the celestial upbore 

SI 


m 

a 


A A 

e„ + ie , 

2 a 


when* e^ 


A 

and c,j 


arc unit vectors and the celestial sphere is treated as an 


almost-complox two-dimensional manifold. The set of vectors forma a complex 
null tetrad from which the metric g ^ is composed; 


’ah 


2 / 


(a"h) 


- 2m m 

(a n I 


•1 *1 »l • »l 

f n « n t' - -m m = -m m‘ -1 

a a a a 


(all other inner products vanish; the parenthesis, ( ), denotes symmetrimtion 
and the overbar is complex conjugate). Then 


I 


lb 

• '.iib" 1 m 


t 


and (he Newman-l*c»roso equation in 


2 - r 

-p - n a - R 


neb 


/ a / h 


( t is tin* shear: 


r t 


a b 

I ...mm 
a I 1 b 


2 2 

For nuil p , this indicates that drt/ds -fl - a - R where s is proper time 
along the geodesic, implying that once n < 0 it is always so. For a static 
metric, tho region wherein 0 < 0 is hounded by a "trapped surface" within 
which all geodesic congruences converge to a singularity of zero volume, 
infinite matter density, and infinite tidal forces 12|. This is a black hole. 

We consider the static spherically symmetric metrics in isotropic form 



, 2 V 


+ e 


(< 


2 9 *> 2 2 2 

nr + H ir niiTflcty ) 


(i) 


where 


0 * 0(H) 


and 


V = i/(H) . 


Wo use tho isotropic form ulnoo mot* ! closed-form solutions for various gravity 
tlioorioH are given tin* isotropic form. Wo trnnslato to a luminosity (retardcd- 
tin».;) coordinate via 


o*du 


c^lt - o'^lM 


to find 



o*^iiu 2 - *o***ih*m + o* l| ii 2 <ifi 2 


( 2 ) 


where 


2 2 2 2 
d«“ » dfl“ + sin fldd~ 


From the preceding discussion we find 


P - (1 + , (3) 

whore Di/YDH . Since there is no rotation, p = 0 and 0 < 0 if 1 + i/Tl < 0, 
where n= 0 at the trapped-surfaco radius lt ( . This is oloarly so if t! ean be 

represented as an inverse power series in H, 


* c S a n ( i/R ) n 


for integer n and non-negative constants a . 

n 

Examining geodesics, we must sec if particles in orbit or incident from 
a great distance arc captured by the hole, because a black hole must consume 
anything that comes too close. For metric (1) the geodesic equations with 
motion in the equatorial plane give 


i 


n 




h 


ami 



E . 


two conHtantH of the motion defining the impact parameter X = K/h . Wo can 
then arrive at the equation 


dR/d* « ±(R/X) IH 2 d 2 ^^ -X 2 + (« 2 X 2 o 2(, R 2 )/h 2 | 12 


( 4 ) 


where ( °,-l for null and timeliko He cause material particles usunlly 

have a more difficult time of it than the photons, wo need only look at e = 0. 
Then, 


dR/d* *(R/X) (rV^"^ -X 2 1 1/2 . (5) 

Orbits are stable down to a critical radius R Round by solving dR/d+ 0. 

c 

The capture impact parameter iB X , found from d/dH(dIt/d*t») ■ 0. For R 

and X this Rivets 
c 

1 + R(*' + *’) I - 0 («) 

c 

and 


X c = R c OKp(2|0(R c ) -<f(R c ))} . 


( 7 ) 
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The* results an 1 given in Table 1. The results for the Iteissner-Nordstrom 
metric in Finstein's theory nro presented first. Listed across the table are the 

2* oj 

theory designations e"' , e“ , R^, and r { , the trapped -surface radii in isotropic 
and schwarzschild coordinates 



f. ,2 

-e dl + 


c 2 
o"dr + 


rdst 


2 


the two capture radii It (isotropic) and r (Schwarzchild) , and the 

c o 

coordinate -independent \ . 

c 

This obviously disproves the previous claims that black holes do not oxist 
in Rosen's theory (5) or are not approachable by a geodesic ( 22 1 in the I.ightman- 
Lee theory (('»). Black holes nro, indcod, found to exist in the Kilrnister-Yang 
theory [7|. In another solution for that theory (23|, 


2<p 24 

o=e 


( 1 - m/R) 2 


and 


p = R/( It - M) , 


which is an impenetrable barrier at It M or r - 0 . The Rrans-I)icke [■!! 
black-hole surface has imaginary radius unless u < 0, which is a viable form 
of the theory [ 24 1 • This disproves all prior claims that llrans-Dicke black holes 
are like those of general relativity, "Schwarzschild" [2. r »|. Clearly, all these 
theories predict static black holes, whether or not their proponents have 
claimed so. 

For black-hole evaporation wo follow the discussions of Davies, Fulling, 
and I'nruh (20). We choose the scalar field in a two-dimensional space time 
(of signature zero) because the quantum field theory is solvable. The metric of 
such a space time may always be written in a conformally flat form 


TABLE 1. RESULTS ON BLACK-HOLE EXISTENCE FOR SEVERAL PRESENTLY’ 
KNOWN VIABLE GRAVITATION THEORIES 




















ds m C(u,v) dudv 


whore 


u - t - fo" ' d|{ , 


v = t + /V^dR , 


and 


C(u,v) » 


Using a geodesic point-splitting renormalization procedure, the proponents find 
a satisfactory form for the stress -energy tensor. There are Bevcral coordinate/ 
regularization-dependent parts of 


ah 


«nh R 


dc 


He 


/AHn t 


B ah /N '* 2 * 


ab 


which represent vacuum polarization contributions, where H ^ is the 

curvature and c is the point-splitting parameter. The first two terms may be 
regarded as hydrodynamic and cosmological constant contributions; 


0 ah = -(12 jt)“ 1 C 1 2 0 2 (C _1/2 )/0x a 8x b 


is tliC vacuum-polarization/scalar particle-creation contrilnition. 

The apjx'arance of radiation depends on the boundary conditions imposed 
on the coordinate system at infinity 127,2s). l or the collapse situation, this 
means that the surface at asymptotically flat infinity must accelerate from the 
surface of a thin shell of matter which provides the mass M in the metric 


10 


C(u,v) (equivalent to tlu* shell accelerating away from that surface as it 
collapses). This implies a coordinate transformation 


u Kfn( A - u) 


where K and A arc constants depending on the mass and the particular 
collapse dynamics. We then use 


O _ 

ds~ C(u, v)dudv 


in the T ^ computation. Ilecause such Ixiund.iry conditions may In- found for 

any time-like 2-surfacc 12(»|, the collapse T ^ can lx; determined by the same 

formula for all the gravity theories previously discussed. We have computed 

2 

ft . and found a constant Hawking flux term of the form ~1/K for all the 

al) 

metrics discussed, it extracts mass from the black hole as particles are 
radiated to infinity 127). Thus, black-hole evaporation is very likely a common 
feature of all gravitation theories. 

This establishes that black holes are a normal, not pathological, feature 
of a viablo gravitation theory. \lso, when gravity, thermodynamics, and 
quantum field theory are properly married (menage a trois?), static black 
holes evaporate, radiating particles. Such a process may even prevent their 
formation [21|. Thus, it is useless to use the existence and cv(i|x>ration of 
black holes as a test to determine tlx- relative viability of gravitation theories. 


III. NONVIABILITY OF THE TWO-TENSOR PLUS SCALAR 
GRAVITY THEORY OF YILMAZ 

A number of scalar-tensor [2!>| and two-tensor |:jo,31| gravitation 
theories are currently viable. In certain of the two -tensor theories, the trace 
of one of the tensors plays the role of a scalar field. Alternatively, one could 
treat tin* scalar field of a scalar -tensor theory as the trace of a second tensor 
and generalize the theory to include the second tensor. Cine such theory [s| 
and its generalizations (0, 10 1 are analyzed in this section. The results of such 


I 
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generalizations are shown to pose severe problems of consistency, uniqueness, 
covariance, and coupling of gravitation'll fields to matter so that gravitational 
effects exist. Thus, unfortunately, we find that this theory anti its generaliza- 
tions are not viablo. 

The 195‘> theory [8| leads to the static metric 


.2 2M/r 2 2. 2. -2M/r.,2 

ds » e (Ir + r d^ ) - o dt 


It is essentially a scalar-tensor theory. The results satisfy the three classical 
tests. Pago and Tupper corrected a Langrangian error for the theory and 
showed that the corrected theory possesses non-unique field equations and that 
the static metric given previously is not the unique spherically symmetric static 
solution. A similar theory of Papajietrou has been shown nonviablo. Ni has 
shown that it also is in violent conflict with a few experiments, including the 
perihelion shift. The subsequent development of the theory to a two-tensor 
theory has failed to correct these problems. 

The second tensor h . generated from the scalar as its trace is brought 

ah 

into the theory [9,10| via a local Lorentz invariance argument and related to 
the stress energy tensor via the equation 


lla 


'ah II a 


4nT 


ab 


where II is the covariant derivative. The metric g ^ is then found from the 
relation 



2 


A kind of Einstein field equation is then given by 



8 ' <T ,b + V 4 ') 


12 


These terms 


where t . is the energy tensor for the gravitational field h . . 
ab an 

are as given in the references which report the development of the theory [9, 10 J . 

A gauge condition 




0 


leads to 



-2 U. h d 

V dla c 


Ih 


- oK 


ib h edle h 


•die 


h |a h |b 


‘ B ;,b h 



The geodesic equations are imposed to describe matter trajectories, as is usual 
in attempts to formulate a gravitation theory embedded in a Reimannlan mani- 
fold. One finds 


du./ds 

j 



k 

u u 
l 


u. = dx, x ds 

J J 


The metric g is related to the Lorentz metric h and second tensor h by 


g = n • exp| 2(h 1 - 2h) | 


The metric g to third order in h is then 



^ ( h *n - tih“h . + 12hh.‘ l |i . - Hh ‘'li . h. h ) + n. + ?(hn. - 2h. t ) 
3 \ ij ij i aj i ah j / ij ij ij' 


+ 2 




Ih. .h + ih.'h . 

ij i aj 
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Questions of uniqueness, self-consistency, and actual predicted effects 
arc examined in this theory. First, tlie field equations involving the Einstein 
tensor differ in third order from the previously cited expansion. In fact, we 
find on substituting it into the field equation that the? field creations do not imply 
the relation butween g , n , and h given previously to third order. Therefore, 
the theory's consistency ends there. 

Second, in looking for useful approximations, attempts have been made 
to iterate the field equation to third and higher order to prove integrability of 
the field equation. For parameter A « 1 , we substitute 




"ij + 


2 3 4 

Ap . + X k, , + A .1 4 C)( X ) 

iJ ‘J iJ 


into the field equations and find 


K 


ij 


n..( 1 ♦ 2h 4- 2h‘ _ ) - 1h. - Hhh.. + Hn .h^h* 1 

ij' ' ij ij ab i j 


up to second order consistently for p , and k . Hut we also find 

»J »j 


J., - ^(h'V “ 6h~h 4 ^hh^h . -24 /h.^h.^dh , 

ij 3 \ ij ij i aj J i j ab/ 


showing that the equations are not integrable to third order. Thus, the equations 
can only be; consistent and integrable to second order in h ^ . However, one 
can go further by rewriting the field equations 




II a 


2t.. 

i] 


f a ) ia*°< h3 > 


0 


E 


a 

ij 




.a . b 

h bl ( i h j ) 



i h j»lh 
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To verify, the gauge conditions must be used frequently. E‘ ... does not 

ij I a 

vanish; therefore, the field equations are themselves inconsistent to second 
order. 

Considering t and using the gauge conditions, we find a large portion 
of t ^ vanishes. Therefore, it is not the energy density of a spin-two field and 

t| * 0 ; consequently, Noether's theory is violated. At this stage it is clear 

that this is at best second order theory. Continuing, we use the h field 
equation to write the geodesics as ,l 


pdu./ds = p|, j k | “j uk = 


k II a 


*J II a 


where p is the matter density. The Christoffel symbol to first order is 


I a 


( 26 V 


= 2 [ 26" ,,h . . - n. .h 1 a - 2h . 1 - lh 

Ij) ij ij 


;«) 


Insertion of this into the geodesic equation, partial integration, and use of the 
gauge conditions consistently leaves 


i J k “/ = ° ' 


and the equation of motion is 


2 . 2 
d x./ds = 0 
i 


Thus, there are no gravitational effects in the theory to all orders. All motion 
is simple inertial motion in the theory; therefore, there are no gravitational 
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effects. We conclude that one must bo very careful jn generalizing a scalar 
theory to a tensor theory. In particular, gauge conditions must lx* kept an 
auxiliary ( not essential) part of the theory. 


IV. HOMOGENEOUS ANISOTROPIC COSMOLOGICAL MODEL 
IN ROSEN'S GRAVITATION THEORY 


Rosen's theory is a two-metric theory with prior geometry (5). We shall 
consider one of the metrics to in? tin? I.orentz metric y n . The field 

' a , n 1 * 

equations in Rosen's theory for g are 


T n fl 'aft 


N 

pi' 



-Kirk T 

pe 


where 


N * 
pe 



and 


k 





The slash (I) indicates the covariant derivative with respect to 
is the stress density tensor. 


y , anil T 
pe pi' 


Since g is a Ricmannian metric, we adopt the Rianchi classification 
pe 

scheme to the homogeneous three spaces [32|. We also wish to express the 
lield equations in terms of the irreducible parts of the stress tensor which 
possess invariant physical significance. Hence, one desires to have equations 
k j 1 i l\ 

for T , T . , T. , and T -—6. T, (the I.atin indices running 1-3). 
oo oi k j 3 j k ' 


10 










If the metric 1 b expressed in terms of a basis tetrad 


B - E A E a , 

ftV M \V 


then the affine connection coefficients become 


OT 


E A,n 


I 






and we have 


R , - T v ♦ r p 

p r I a fia vn 




Then Rosen's tensor N may be constructed in terms of the affine connection. 


■ ? r,,x (Ci« * r li« ) -* (C * r * )( 


r 1 ' + r x 

Xn 


vn 


)■ 


and the field equations written in terms of the connection coefficients implied by 
a particular set of symmetries chosen for the spacetime. The Rianchi classi- 
fications may then bo used. 

A suitable metric would bo Riven by the line element 


. 2 2S2 . 2 2o 2ft _i J . A . B 

ds = e dt - c e . . E^ E^ dx dx , 


where n, o , and ft^ are functions of the time. The metric is similar to a 
form widely used now in Einstein's theory (32|. 
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I 






The Hi unchi case 


P i , i 
e a ■ 5 a 


will now lx* considered. The dynamical field equations are then 


-H, r< 


V^V 4 (v3 4 M ( ) 




where r\. are the traco-froo stresses. For a diagonal /<.. , the fluid shear is 


"u ■ % • 


The shear field equation 




where X = viscosity is then 




-HttAiTj. - Ho 0 . . 


and for small shear (o' - « 1) we find 




[const. ) o 



-H?rX(t-t ) -8(o-r> ) 

o 


Anisotropy damping is much more rapid in Roaen's theory than in Einstein's 
theory. 


For a dust source?, the hydrodynamic conservation law wives the energy 
density as 


-3ev 

P ■ P D ° t 

as in Einstein’s theory. Then, we have 


7:0 + 


Gri 2 


„ 2 -Uia 
2 a o 
o 


- 8 m) e 
o 


-:v> 




and from the 



equation 


• • 

rv 


may be eliminated to find a Friedmann-like 


equation. The Itinnchi classifications are being investigated in this context. 
The most interesting asjxjct of this is the isotropization question: VK ill Itosen's 
theory provide a bettor key to the isotropy of the universe? 


V. INHOMOGENEOUS COSMOLOG I CAL MODELS 
IN EINSTEIN'S THEORY 

A. Tilted Self-Similar Cosmologies 


Self-similar cosmologies are models that admit a group of similarity 
transformations; i.e., an invariance under changes of length scale [J5|. There 
is defined a homothotic vector by the relation 


V 
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where X is a scalar ami V is a covariant derivative. Tilted models are 

M 

those in which the fluid flow voctor is tilted away from the time -like vector 
normal to the space-like hy|xjrsurfaces rj by an hyperbolic angle (i [33]. 

We combine the two formalisms in search of interesting models. The particular 
question is; Arc there tilted universes with submanifolds, invariant under the 
Hianchi 1 group (translations), with those manifolds submamfolded to ones 
invariant under a similar group and with the fluid flow shear free? Kardley has 
mapped out the self-similar formalism 115J, and King and Kills describe the 
methods needed for tilted models (33). 

The equations of Kardley were modified to include a tilted flow vector and 
specialized to the case with lUanchi I submanifold. The undefined quantities 
are those given in Kardley 's paper. The metric is 


”"" b ] 


where z labels invariant submanifolds and di, - b a (b is a voctor related 

a ' a , 

to the structure functions of the similarity groiji H ). The generate a Cl 

u Z 

subgroup characterized by the Hianchi classification, and 


B ab< 7 > 


2rv 2fl 
e e . 
ab 


a » o(z) , 


'’ab * "nb< 7 > • 
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For Ilianchi 1 submanifolds invariant under (the field with shear-free 
flow), the field equations are: 


-U 2 *b 2 - T° o' 2 * 
o 


. o -till 

•V • V 


n ■ *- ... - rpk -2li 

-tin - 9r> + 3b » I , e 

k 


2b b - - 6 b 2 » ?r o 
it 3 if if 


where 


dx/dz 


Reducing T to its irreducible parts, a perfect fluid source, 
\iv 


T = ( p + p) u u + pg , 

lit* H V (A I' 


is obtained. With further specialization, the field equations and fluid conserva- 
tion laws (shear free — p 0) are: 

.2 -2rv i T 2 *1 -2i I ORIGINAL 
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Km ♦ p) Vo 1 




-2e b. a 
k 


v 2 *' 


*!*• • 2 n "2a , . i 

-Go - 9rv ♦ 3o bjb 



f 


5 6 „ b \ 

-2i 

-2a e 

e - — ■ 

m m 

a k 1 . k 

(n + pluu - — o u. u 
\P M j u f 3 k 

. 


. _ 


b ♦ 3 (p + p) a ■ 0 


(„*p)u o - V p |o 
and 

(p + p) &k - -V°P| 0 


The equation of state is p yp an<l Uie equation for flow normalisation is 

u i/ 1 ■ -1 . 

M 


For the hydrcKlvnamic equations, 


p - P Q expl-3( 1 + y) nr 1 


and 


u = u oxp( -3-)fY j , 
a a 


ami for the trace -free stress equations 


. r~o /iA . . 1 1/2 -4 -3(l + 3v)rt/2 

b i " M 1 +V) P</ 2 0 • 

Tho Friedman-like equation for a in solved after the substitution o « fn a 
to Rivo 


u (h ♦ i) 'i 
(3v+l) 


a 


a 

o 


z 

X I o 1 dz . 

O •' 

z 

o 


Tho density, p, c:m lw arbitrarily specified as a function of the space coordinates 

at time z , p p (x) . Since these are self-similar models | 31 1, z= t n x n ' , 
o o o 

where n and m are real numbers and x is any space coordinate. The function 
1 / is fixed by the condition <&,' | ^Tz> ^ 0 (15). 

Since 

dip m b o 1 = b dx*' , 
a a 


suppose we choose 


b 1 0 ; tlicn wb do, indeed, have the solution 

•f 


3 a 


(3y+ l)/2 


-x 


/ ( 3y+ 1) 


x e 
o 


3 n m 
1 *2 


B. Universes with Nonsynchronous Time 

Now, let the metric lie of the form 




u z 
a 


a^2 


v‘ z) 


t 


ORIGINAL PAGK L 
OF POOH QUALl 1 V 


23 


where the forms 


z a » A(t) B a (t) e| dx* 


■ e 


« J 

a . a 
aj 




and where 


“a * “a*'-** * 

is a function of time and the space coordinates. 
In this system, the velocity form is 


a o 
u = -z 


u 

a 



and the "time" coordinate is 


. o .. a 
dx = dt - u 7 . 

a 


This makes a very convenient system for computations; i.e., the projection 
operator is simply 


,v . V , O - P 

h =6+66 . 

\i n n o 

The convective derivative? , u'V , is given as V and, in the case of scalars, 
o a ° 

is given as 0/9x . The kincmntical gradients of the fluid velocity arc also 


21 




easy to compute Bince tho ! J (purely spatial) rotation coefficients Ho not 

o 

appear in tho computations when u^ « - 6 ^ . 

In terms of acceleration a., expansion S , rotation W , reference 

i mr mr 

frame spin t , and proup structure functions cj^ , we find 


.i 

oo 


r__ ■ a. , 


r 


m 

ro 


8 + W , 

mr mr 


r 


r 

om 


rr 


-t 

mr 




I 


1 

jk 



l ik U j 


* 



for the rotation coefficients. In terms of the metric functions, these kinemntical 
quantities are: 


S 06 + a + u. , . + u u 

mr mr mr (mlr) (m r» 


W 

mr 


u, , . + u r u , + ut , u , + u, C k 
(mlr| [ m r| jj[m r] k mr 
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and 


u i * V) 


where 


mr 


, ’k<m n r> k 


Acn^rlk 


t 

mr 


v -1 
B. . B 
k| m r 


1 k 


(</ ^kf m^ C \]k 


and 


t 

mr 


n, , n"! . 

k[ m r ] k 


* e() k[m (e \]k 


ri 

If we write the shear tensor > as 

-'mn 



S 


mn 


1 

3 


0 

mn 


0 


and the expansion scalar o as 


0 



% 


then the field equations for submanifolds of Bianchi I are: 



IV V + I r 

2 -'mn i-mn 2 


=s 


p 
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( l<,8 km> * 4, m ) “ 


kSm t /m U m 




and 


£* +M ^nk* IMI nk* r 


nk 


- ^ 6 .r 

3 nk 


nk 


whore y ■ 8 y/0 x° ; the spatial curvature r is caused by the t u term in the 

1J K 

rotation coefficients; q, and it are the momentum-density flux and trace-free 

I UK 

stresses. The hydrodynamic equations are 


p + (p + p) 0 


and 


(p + p) u = -h” p. 

' m m In 


The reader may have realized by now that the purpose of this discussion 

is to demonstrate the tri virility of finding inhomogeneous solutions. Therefore, 

we immediately specialize to the seemingly most unpleasant case: Perfect fluid 

( it , = q. 0), diagon;d p , and dust (p^- 0). Then, the acceleration 
nk / mn 

equations read 


a 

m 


u 

m 


+ u. 




0 


which solve to 
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At "'O' ~P 

o ni 

u « u (x) e 
m m 


o 

where u (x) is a function of the space coordinates. Now, a and a are 
m o o m 

functions of x , o(x ), and /j (x ) . Let a =* Ina and the Friedman-like 

equation is 


da/dx 


1/2 


P ° 

4 + - 2 - 

. 3 „ C 

La 2a 


o m 

a - u . e 
mlm 


Since a = a(x°) , we are free to find a solution with u° . =0. Then 

' mlm 


m 


V 0 
m 


where 


* - 2>„° n (I) 


and a^ is a constant determined by the boundary conditions inserted into the 
usual Fourier-Bossel integrals 


a = ff(x) 0 (x)dx 


and O n (\) is any set of orthogonal functions consistent with a chosen symmetry 
for the flat -space Laplacian implied by 


_ m o _2 
V u 0 = 0 . 

m 


2H 




An implicit integral is easily found for a(x°). The explicit dependence of 

|| 

n lx ) is found from the shear equations, 
m ' 

C. Extended Locally Rotational ly Symmetric 
Class 1 1 1 Cosmologies 

Wainwright and Szafron (17(35, 30] have extended the locally rotationolly 
symmetric (I, US) universes |16| to inhomogeneous models of the Sckeros |37J 
type. They suggest a general adaptation of Uiese techniques to study large 
classes of solutions to Einstein's equations with a dust source (17,30). It a the r 
than concentrate on a hack-typo duplication of Class 11 results for Class 111 LRS 
(although wo shall do so eventually) , refer to Paragraphs V.A and V. H. 

Finding outrageous, inhomogeneous solutions is not a very formidable 
task, as previously revealed in Paragraphs V'. A and V. II. The arbitrariness 
allows any density's spatial de|x;ndencc to presume spatial hypcrsurfnccs. 
Therefore, inhomogeneous and single solutions are neither very challenging 
nor very informative. 

Szafron, Wainwright, and Eardley have suggested a solution to tin- pre- 
viously mentioned problem. ItaUier than groping solution-by-solution through 
Einstein's equations, group structures should instead lie used to give an elegant 
classification to families of solutions and relate those families to each other 
mainly as each other's subsets in the group parameter space. 

It seems that these inhomogeneou- models are of little use individually. 
Also, the solutions seem too easy to obtain. An "easy" way out would lx* to 
compute the observational calculations in each case. Them, -me would test for 
• • <> ■ ’ : i wed Inhomogeneoui coi mological modi Lb in comparison with observations. 
However, tins procedure scorns tedious and unrewarding. 

It seems better to study the structure of various generalizations from the 
homogeneous cases to the inhomogeneous cases. < >nc then can treat the problem 
as die consideration of successively larger classes of solutions. The main con- 
sideration is the alteration of the group structure of the submanifold isometric 
under the lower dimensional group. One wants to know what physics there is 
corresponding to the choice of the higher dimensional group of which the addi- 
tional free parameters will allow characterization of the enveloping inhomogeneous 
models. The only available method would lx* the* reverse of the Inonu-Wigncr 
contraction ( 14| . 
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In outline, if we have a group of generators 
algebra 


I with a particular 


'V v 


c** 1 I 

lk 


and the group is not compact, it may be possible to find a largcr-dimension group 
to which the I group is a continuous subgroup. We seek a group of generators 

J. which are found from I by a singular transformation 
lv lv 


lv 


I + cV I 


lv 


HV 1 


and 



where 


V 

HV 


is a matrix representation of the group and as 


c 


0, 


J 


lv 



Then the larger algebra is 


lk -1 2k 

[J, ) « C“ J , . + f C, . J_. = (I. I, | 

lv V li', 1/i lk ly, 1 p 2k lv , lji 


+ e 


v / 6 / +6 , V , + eV , v * |[ 1 1 • 1 1 1 • 

vv jip vv mm vv MM I * v *M 


■) 


30 


and as e — 0 wc recover the subalgcbra. There is as yet no elegant way of 
performing this reverse process. Possibly, examination of the group parameters 
rather *han guessing at the representation V _ of the higher dimensional group 
might work. For this method 


M v 


« b +eV b/^ 

li' \v ip 1 


and 



These parameters are directly related to certain relative tensors which define 

the structure constants cj. of the algebra ( 32 ) . At the same time they 

jk 

directly specify the vector fields in the manifold which generate the algebra. 


Spero and Baierlain | If) ] have suggested a variational approach which 
may lead to techniques for solving this problem. They treat the inhomogeneous 
metric as possessing an approximate symmetry. However, the following 
question should be considered when applying this solution: Under what conditions 
is this a "subsymmotry" of a higher-dimensional symmetry characterizing tlx; 
inhomogeneous model ? 


Tlic Class 111 I.HS models offer the most challenging testing ground and 
shall be considered in a future publication. 


VI. HOMOGENEOUS COSMOLOGICAL MODEL IN VANG'S 
THEORY OF GRAVITATION 


Yang (7| has presented a gauge gravitation theory which is a rederivation 
of Kilmister 138). In this section, a cosmological solution is investigated. 

A metric is chosen in the form 



,2 , . 2 2 , 2 . 2 . 2,2 

dt - Adx - Ax dy - Ax sin ydz 


and implemented in the MACSYMA symbolic manipulation system as the matrix 


31 


-A 

8 

8 

8 


8 

-AX 5 

8 

8 


-AX 2 SIN 2 (Y) 8 


where A is a function of t and x . 

The Kilmlster-Yang equations R j(j n ^ j * 0 then have the following 
surviving four components ( OA/0t ■ A^ # 0 A/OX ■ A^) 


[ 


(D5) | 2A 2 A t A xx X 2 - 2A 2 A TX A X X 2 ♦ 2A 4 A TTT X 2 - 2A 3 A T A TT X 2 

+ 6A 2 A T A X X - 4A‘ , A tx X 


4 3 3 32334 35 3 

4A W - 14A W + 9A A x x +4A a tt a x x - 4A a ttx x 

4 2 3 2 2 4 

+ 8A A^X - 16A A^X - 8A A y X , 


6A 3a t A xx x2 “ 9A 2 A t A x 2 + 8A 3 A TX A X X 2 + ^’’a^X 2 - 4A 4 A T A TT X 2 
- 4a4a txx x2 ♦ 1«A 3 A t A x X - 4A 4 A tx X, 2A*A ttx X 2 - Za’a^A^ 2 


2A CX - 4AA xx X + 3A x X + GA A^X - HAA^ 
2A 3 X 


32 


T 


Assuming C is the constant scalar curvature, A(t,y) is separable, and 
A * I* (X) Q (t), the equations take the form 


(Dd) 


2 3/ 
2P Q X II 


p ,qq ttt x " ,, ,q t q tt x * pp xx q t x " p x 


o x + pp g \ 
T X V T I 


•> 5 / *> 

*~g x 1 4p 


J W 2 - ""’xV* ' «*xV * HpL,|> XX X • ,6PI X X - 8p2p ] 


) 


V* ^ll' W r| | X - 4r : '<|, r Q r| .X 


• 2PP xx q t x - p x Q ~ 


V * "’Vr) 


. » . 


«P 3 Q X + 2CP 3 ()X - 4PP X ♦ 31>^x - 8PP 

1 1 A A A ^ 

2 p 3 gx 


The second of the equations involves only the function 1* (X) 
-1 


A solution is 


I'pon sui.stitution, the following eouations for Q survive 


(w TTT - <V'tt> 


1H 


4Q‘‘(QQ ttt - Q t Q tt ) 3 Q tt +CQ 

* s ' 


X 
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Q 


each component separately equal to zero. 

The equation involving the curvature scalar is 



eg 

3 
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t'pon differentiation, we find 


Q 


TTT 





which identically satisfies the remaining field equations. Thus, we find 


Q = Q o exp l */"-c" 73 t | 


The metric in then 



dt“ - Q 

o 


\' -c /:i t 

o 


-4 2 2 2 

X (dX +X dsi ) 


which can be easily transformed to 



dt" - cxp| \ f -C /3 t| Idr" + r *'«! S2 “ J 


A similar metric has been found as a cosmological solution in a modified version 
of Yilmaz’ theory [39|. It is also similar to the metric in the steady-state 
theory |3|. Therefore, we refer to the solution as the Yang-Yilmaz universe. 

The value of the scalar curvature C determines the large-scale behavior 
of the solution. For C > 0 the universe is an oscillating one, and for C < 0 it 
expands forever. This is a solution which is not an Einstein space. 

The observations in the Yang-Yilmaz universe for C < 0 would be very 
similar to those in the steady-state universe |3|. For C > 0 we have a curious 
mixture of flat spatial sections with a closed universe. For real solutions we can 
arrange so that 


Q = Q cos Kt 
o 


34 


whore 
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K - >4 -Cl 3 . 


Then the universe oscillates in volume with the evolutionary |ieriod 


2*/K - 2tr < 3 / \ r -C . 


Present observational evidence seems to indicate an open universe (the C > 0 
case) (40 1 . Nevertheless, observations would seem to lie inconsistent with 
either choice of C; and, of course, for C 0 we have only Minkowski space. 

VII. TOWARD FINDING THE ROTATING SOLUTIONS 
FOR GRAVITATION THEORIES 


The gyroscope experiment (41) was proposed to test whether an induc- 

t op field (caused by rotating sources) exists in gravitation. This experiment 

is considered the "Faraday" experiment of relativistic gravitation theories. In 

the wcak-field limit, this tost of the existence of such fields is essential. The 

existence of the binary |>ulsar allows extension of the gyroscojje experiment's 

results to strong gravity fields; therefore, it may Is- of use in testing gravitation 

thoories. To do tins the fields exterior to rotating stars must lie found for each 

viable gravitation theory. This may lie |H»sed as the problem of determining the 

axially symmetric, stationary asymptotically flat solutions for empty space for 

these gravity theories. Since most of these theories possess a Iticmannion 

geometry, one possible technique would involve use of the Newman -Penrose 

identities in finding the solutions in which the gravitational field’s directions of 

propagation are expanding and twisting. They may not Ic shear free, and, 

therefore, then; may 1* no analogue of the Kerr |3| solution in Einstein's 

theory. An alternative and more meaningful approach might involve the use of 

the "Killing" vectors of the static spherically symmetric solution to extend to 

the "Killing" vectors of the stationary axially symmetric case as has bee n done 

in Kinstoin’s theory |42|, belying on that, one may make use of the tetrad 

vectors E from which a metric g , is composed, 
ap ab 
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I 





and try to extend 

Tho solutions for tlto nf -it ««• spherically symmetric ease in tlio viable 
gravitation thoorio* aro in tho form 


•» 

flu' 


- o 



'i 

i 


oj, o 

o' d|{ + o 


2d 2 2 ° ' 2 2 2 

lt~d«* ■» o’ It Min*’fldi|» 


whore 


0 » 0(1*) 


and 


0 - </(«) • 


The appropriate basin of differential forms u f* are 


( . 2 n v \ o 0 

ds » n u cr ) u’ ■ e dt , 
M" / 


1 i / 

u> = edit , 


2 J 
w = o Rdfl , 


and 


3 i/ 

us c Ksinfldv 


3fi 


I 


Wo can immediately extend these forms to tiie stur onary axially symmetric case 
by writing I 1 3 1 


u° » e 0 |dt + f ( 0) d'Pl , 


ui • c % dlt , 


2 il 

u • o ltd n , 


and 


3 </ 

u> - o Rsiivid? 


I ho contravariant basis dual to these forms is 


(< a/ 1 , e > - 6 11 ] , c 
\ v v J o 


ej - e^n/on , 


0.^=0 It B/l)f) , 


”0 . 

e » /at , 


and 


c |c It sin# | *0/8 - *(») | e It Hino | 1 e"°o/f)t . 

•I 


The composition of tlie metric is 


K e‘e If Ef dx 1 ) 

V M »'* \ i / 
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anrl 


J lv 




J 1 O/Dx' 


From this wo have 


o it) 

E * c 9 , 
o 


and 


Ky * o 0 f(n) . 


r 1 J' 

Ej * O ( 


2 1/ 

K., = o It 


3 »/ 

E, } = e Hsinfl 


for the (invariant basis, and 








.1 = e 


j: 


-<i‘ 

e , 


2 -i/ 

= c /K , 
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■ (e^R ainfll 

and 

J 3 * -f(fl)“/lc H sinfl | . 

o 

It is this vector decomposition that should be inserted directly into the field 
equations of the particular gravity theory and the solutions found for the functions 

E** and .1* . The Newninn-Penrose identities can then bo used to investigate 
i M 

the properties of the solutions. In view of the fact that for all the alternative 
theories of gravity H * 0 , it may not lie jxissible to find the solution with 

shear-free expanding and twisting rays; hence, there is no analogue to the 
Kerr I'll) or Newman-Tamburino-l’nti (45J solutions in Einstein’s theories. 
These solutions are presently l»eing worked out, and it is expected that they 
will be obtained systematically and with facility. 

The equivalent in matter-filled spaces would be shearing universes with 
twisting rays. These solutions also will lx; investigated using the algorithm 
previously mentioned. 
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